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SIMULATION  OF  UlOr.RKATE?  DIVERSITY  ARRAY 


In  earlier  reports  [1,2]  a  three-dimensional  underwater  array  './Inch  simul¬ 
taneously  forms  and  combines  multiple  beans  in  elevation,  v;as  described  and 
analysed.  -•  The  objective  of  the  system  js  t  o  approach  a  condition  in  vMolt  the 
multiple  ray  arrivals  from  a  distant  source  are  separately  received  and  coherently 
combined.  Because  the  analytical  forms  Riving  the  final  output  sign  a1 -tc~noi sc 
ratio  (DVD)  are  too  involved  fur  direct  ct riput  at  ioi^  s  i-vnlat  Ion  experiments  were 


carried  out.  Results  of  thi 


in  { 1 ,  section  3,  pp.  26-35i 


ir.  work  are  presented  fro lavl,  for  thof  system  descri 
b  v.'hi  eh  is  based  on  maximal  ratio  'combining  (MKC)  of 


scribed 


the  nu: It  t:>lu  ravs  formed.  At  the  same  t  ime  certain  other  cases-  v.we  simulated 


[1]  Fred  l'aber ,  "Research  dr;  I- i  :  t  r  United  I’e.  lr  r.vat  •  t  Acoustic  Airiyr.,"  Vi'RC 
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in  order  to  see  what  magnitude  of  improvement  is  obtained  using  this  method  over 
simpler  ones.  Given  below  are  results  for  (1)  a  three-dimensional  array  with 
multiple  beams  (branches)  using  only  coherent  phasing  of  branches  (known  as 
equal  gain  combining  (EGC) ,  (2)  a  three-dimensional  array  with  multiple  beams 
using  selection  of  the  maximum  amplitude  branch  (known  as  selection  combining 
(SC)),  (3)  a  three-dimensional  array  with  a  single  fixed  focus  beam,  and  (4)  a 
two-dimensional  array. ^_The  last  one  represents  the  original  concept  explored 
earlier  and  reported  fb'i3,4,5]. 

The  simulation  results  can  be  briefly  summarized  as  follows.  For  the  con¬ 
ditions  chosen,  an  improvement  of  approximately  4  to  1  is  obtained  using  the 
maximal  ratio  combining  technique  with  the  three-dimensional  array  over  anything 
else  that  was  done.  An  interesting,  though  expected  result,  was  that  the  var¬ 
iance  relative  to  the  mean  of  the  output  CNR  using  this  best  technique  was  much 
less  than  for  the  two-dimensional  array.  The  computational  model  assumed  ten 
independent  ray  arrivals  randomly  spread  over  ±10°  jn  the  vertical.  In  the  two- 
dimensional  array  these  combine  noncoheroutly  resulting  in  a  nearly  Rayleigh 
fluctuation  of  amplitude.  In  the  three-dimensional  array  the  ten  rays  are  to 
a  large  extent  resolved  in  the  separate  branches.  The  diversity  selection  or 
coherent  combining  then  results  in  a  substantially  smaller  fluctuation. 

We  now  describe  the  physical  and  statistical  arrangements  assumed.  Fig¬ 
ure  2.1  suggests  the  deployment  of  elements  in  a  three-dimensional  space.  Element 

positions  (X  ,Y  ,Z  >,  rs  =  1,2  ...  N  were  assumed  independent  random  vectors,  the 
n  n  u 

number  of  elements  X  being  31  for  this  computation.  The  horizontal  coordinates 
(X  ,Y  )  were  assumed  independent  normally  distributed  random  variables  with  zero 
mean  and  .standard  deviation  of  50  wavelengths.  This  value  of  standard  deviation 
implies  that  at  100  Hz  where  the  wavelength  is  about.  15  meters,  about  68%  of  the 
array  elements  will  be  concentrated  in  a  range  of  1'  750  meters  around  the  center 
of  the  array.  Because  the  array  main  beam  was  focused  to  look  in  the  Y-Z  plane 
only,  the  random  variables  X  were  not  involved  in  the  computation.  The  vertical 
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FIGURE  2.1  PHYSICAL  ARRANGEMENT  ASSUMED  IN'  SIMULATION. 


coordinate  2  was  assumed  unifonr.lv  distributed  with  trie  an  zero  and  range  of 
n 

50  wavelengths,  also.  This  implies  a  depth  range  of  t  375  meters  around  the 
array  center. 

Ray  arrivals  were  assumed  to  be*  in  the  Y-Z  plane  corresponding  to  the  azi¬ 
muthal  angle  of  focus  of  the  array,  but  the  rays;  were  assumed  dispersed  in 
vertical  angle.  Tan  rays  were  as  surra  d  arriving,  all  of  equal  magnitude  and 
independent  random  electrical  phases  :<  ,  m  -  1,2  ...  10,  each  uniformly  dis- 

IU 

tributed  in  (0.  7i).  The  arrival  angles  0  ,  measured  from  the  vertical  were 
assumed  independent  and  uniformly  distributed  over  (80p,  100°),  or  _  10°  rela¬ 
tive  to  the  horizontal. 

Ten  sets  of  random  pairs  of  numbers  (y  ,?.  )  n  =  1,2  ...  31  to  represent 

ten  possible  random  element  positions  were  chosen.  For  each  of  these  positions, 

five  sets  of  angle  pairs  (3  ,«  ),  m  r-  1,2  ...  10  were  randomly  selected. 

mm 

Assuming  the  combining  scheme  of  [1,  sect  ion  3,  Figeic  3.2),  the  output  of 
ea~h  branch  prior  to  weighting  and  combining  is  given  by 


M  N  i[k(x  sin 0  cor,-''  -Iv  sin:‘  +;-■  cosO  )+£  .-to  ] 

..  r  v  ,  m  m  m  'u  m  m  n  m  nx  m 

vi =  \  i  v 

m-1  n-1 


-  AJc',,i,  i  *  •  1,2  ...  I 


>  (IT 


where  6  )  are  the  azimuth  and  elevation  angles  of  the  mth  arriving  ray,  and 

(Bm,  a^)  are  the  magnitude  and  electrical  phase  at  the  center  of  coordinates  of 
that  ray.  is  the  injected  electrical  angle  to  focus  the  array  to  azimuth  and 

elevation  angles  ($>g,  0g^)  as  shown  in  [1]  Figure  3.2,  page  31,  and  is  given  by 

$  ,  *=  -k[x  sin0  .  cosd  +y  sin9  . s J n d»  +z  cosO  ,]  (2) 

ni  n  si  s  si  ’s  n  si  v  ' 


The  subscript  i  identifies  the  beam  or  branch  number  and  I  is  the  number  of  branches 

used;  here  we  will  use  I  =  17  branches  for  reasons  to  be  stated  below.  (1)  and  (2) 

are  written  to  include  rays  arriving  from  any  azimuth  and  elevation  (a  ,  e  )  and 

focused  at  any  azimuth  and  elevation  (4>g,  6 s^) .  We  now  specialize  these  expressions 

to  the  case  of  the  array  focused  in  the  Y-Z  plane  and  ray  arrivals  in  the  Y-Z  plane. 

Thus  <t>  =  <p  =  90°  and  (1)  and  (2)  reduce  to 

m  s 


M 

h-  I 


N 

1 

m=l  n=l 


_  j  [ky  smO  +z  eosO  ) H  a  ] 
B  eJ  Jn  m  r.  m  Yni  in 
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(3) 


and 


na 


=  -k[v  sin? 

n  f  i 


:  tosO  .] 

n  si 
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The  final  signal  output  using  KRC  is  given  by 

I  ~  I 

S  =  l  |V.|2  -  l  A? 


I-?.  1 


(5) 


The  mean  square  value  of  the  noise  at  the  output  of  the  system  is  given  by 


N2  - 


N 

5: 

n=l 


2 

n  > 
n 


I  I 

I  1 

i=l  j-1 


A,A.cos(<>  ,-<p 
i  j  ni 


(6) 


where  the  (A.,  'l> .)  are  the  measured  signal  amplitude  and  phase  on  the  ith  branch 
x  1  2 

given  in  (1),  the  £  ^  are  given  by  (A),  and  <n^>  is  the  equivalent  mean  square 
value  of  the  noise  input  generated  by  the  nth  sensor.  In  the  computation  <n^> 
was  set  equal  to  1.  Finally  the  computed  output  SNR  is  given  by 
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S2 

SNR  =  ^  (7) 

N 

As  discussed  in  [2]  the  angles  0  . ,  if  set  at  values  separated  by  , 

si  n 

n  an  integer,  will  result  in  uncorrelated  noise  voltages  at  the  branch  output  pro¬ 
vided  the  elements  are  uniformly  distributed  in  depth  over  the  range  (-h,  h) . 

For  the  geometric  conditions  used  here  it  was  determined  that  1.15°  spacing  between 
the  0  ^  would  accomplish  this  when  h  =  25X.  We  have  accordingly  assumed  17  beams 
symmetrically  placed  around  0  =  90°  at  1.15°  intervals.  Separation  is  about  a 
beamwidth  in  this  case  and  the  total  coverage  in  vertical  angle  is  close  to  ±10° 
with  respect  to  the  horizontal. 

Results  of  the  simulation  using  the  maximal  ratio  combining  technique  are 
shown  in  Table  1.  Here  we  show  the  average  and  standard  deviation  of  the  SNR  given 
by  (7)  over  the  five  sets  of  paired  values  (e  ,  am)  m  -  1,  2,  ...  10,  for  each  of 
ten  sets  of  position  samples  (yn>  zn)  n  =  1,  2,  ...  31.  These  statistics  are  de¬ 
noted  <SNR>„.  and  a(SNR)_..  Then  the  50  results  of  SNR  (five  sets  over  (0  ,  a  ) 

Q  <p  Q(p  m  m 

times  the  ten  sets  over  (y  ,  z^))  were  treated  ns  a  sample  of  size  50  and  the 

overall  mean  and  standard  deviation  denoted  <SNR>  and  a  respectively  were  determined. 

•  • 

These  were  found  to  be  <SNR^>  =  480  and  a  =  182.  We  point  out  that  if  a  single  ray 
were  assumed  to  impinge  on  a  single  element  the  output  SNR  would  be  unity  and  the 
variance  would  be  zero. 

To  determine  the  effect  of  a  frequency  change  on  the  output  we  have  assumed 
two  different  situations.  In  the  first'  we  assumed  the  angular  separation  between 
beams  held  at  1.15°:  that  is,  0^  was  held  fixed  at  90°  t  n(1.15°),  n  =  0,  1,  2, 

...8.  The  frequency  was  then  changed  by  factors  1/2  and  2.  These  latter  changes 
were  accomolishcd  bv  simaly  changing  the  valuej  of  y  and  z  used  in  the  previous 
calculation  by  the  reciprocal  of  these  same  factors.  The  vertical  beamwidths  be¬ 
come  narrower  at  the  higher  frequency  and  wider  at  the  lower  frequency  but  the 
angular  spacing  between  beams  remains  unchanged.  Thus  the  beams  are  not  optimally 
spaced  resulting  in  either  uncorrelated  noise  in  the  several  branches,  or  in  non¬ 
total  coverage  of  the  vertical  range  within  which  incoming  rays  are  expected.  Result! 
of  these  calcuiations  are  also  shown  in  Table  1  revealing  a  decrease  in  the  over¬ 
all  average  SNR  as  one  might  expect.  These  results  essentially  show  the  sensitivity 
of  the  scheme  to  incorrect  placement  of  the  vertical  beams.  As  we  see,  the  effects 
arc  not  overly  serious,  the  mean  output  at  the  1/2  and  2  times  frequency  points 
being  within  about  80%  of  the  menu  at  the  design  frequency. 
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In  the  second  simulation  of  frequency  effect  the  array  focusing  was  held 
fixed  for  a  100  Hz  sinusoid  by  fixing  the  angles  <{>^.  The  applied  frequency  was 
then  altered  to  99,  99.5  and  100.5  Hz.  As  a  rule  of  thumb  the  bandwidth  of  an 
array  with  fixed  phase  shift  focusing  is  the  inverse  of  the  time  required  for 
the  wave  to  traverse  the  array.  In  this  case  it  would  imply  a  bandwidth  of  the 
order  of  1  Hz,  or  about  1%  of  the  center  frequency.  The  results  of  the  simula¬ 
tion  are  shown  in  Table  2.  The  overall  mean  SNR  is  observed  to  have  fallen  by 
about  3  dB  at  frequencies  100  ±0.5  Hz  from  what  it  was  at  100  Hz,  thus  confirming 
the  rule  of  thumb  on  bandwidth.  The  overall  mean  SNR  at  59  Hz  has  fallen  further, 

the  level  appearing  to  be  about  what  one  gets  when  one  steers  the  azimuthal  focus 

& 

away  from  the  source,  illuminating  the  sidelobes. 

We  point  out  that  the  array  properties  observed  here  are  all  normalized  to 
wavelength  so  that  at  higher  frequencies,  with  the  actual  array  size  l'educed  but 
with  array  size  in  wavelengths  held  constant,  the  bandwidth  would  remain  at  about 
1%.  At  10  kHz  we  expect  a  100  Hz  bandwidth,  a  value  adequate  for  operating  a 
teletype  communication  link.  Furthermore,  at  this  center  frequency  the  array  hor¬ 
izontal  dimension  measured  between  ±o  points  is  .15  meters,  a  dimension  one  might 
envision  for  an  array  suspended  from  a  surface  ship  or  deployed  around  a  submerged 
submarine.  Tie  array  could  therefore  he  useful,  for  underwater  data  communication. 
Furthermore,  our  work  was  based  on  focusing  by  fixed  phasing  of  elements  of  :he 
array.  It  is  the  fixed  phasing  which  limits  the  array  bandwidth.  By  using  con¬ 
trolled  time  delay  networks  at  each  element  involved,  broader  bandvidths  are  achiev¬ 
able  suggesting  the  possibility  of  higher  speed  data  communication,  or  lower  center 
frequency  with  higher  bandwidth. 

The  multibcnm  three-dimensional  array  was  compared  to  four  other  arrangements 
as  follows: 

1.  The  maximal  ratio  combiner  weighting  circuits  which  multiply  each  branch  out- 
*  _  i  m 

put  by  -  A^e  J  i  where  \\  is  given  by  (1)  arc  replaced  by  constant  amplitude 
phase  shifters  e  ''’i,  thus  (5)  and  (6)  become. 

I 

S  “  }  A.  (8) 

i=l  1 


The  r.idelobe  properties  arc  discussed  further  below. 
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Here  only  phase  tracking  is  needed  but  the  mean  SNR  will  not  be  as  good  as  for 
maximal  ratio  combining.  The  term  "equal  gain  combining  (EGC) "  is  used  in  diver¬ 
sity  communication  for  this  arrangement. 

2.  The  maximal  ratio  combiner  was  replaced  by  a  "selection  combiner  (SC)";  that 
is,  one  which  simply  selects  the  output  with  the  maximum  SNR.  This  technique  is 
also  a  standard  scheme  in  communication  diversity  systems.  Whereas  the  maximal 
ratio  combiner  gives  an  output  SNR  which  is  the  sum  of  branch  SNRs,  this  scheme  . 
produces  only  the  maximum  of  the  branch  SNRs.  There  is,  however,  no  need  for 
phase  tracking,  greatly  simplifying  the  processing.  It  is  possible  that 

in  applications  such  as  the  underwater  case  the  problem  is  as  much  variation 
of  vertical  arrival  angle  as  it  is  multipath.  In  effect  then,  the  array  would 
follow  the  variation  in  angle  of  arrival  of  the  maximum  amplitude  ray. 

3.  Output  was  taken  from  one  branch  of  the  multi-beam  array,  the  one  which  fo¬ 
cuses  horizontally  (r  .  =  90°).  The  mirt-ose  of  this  calculation  is  to  see  what 

'  si  • 

effect  is  obtained  when  the  three-dimensional  array  is  operated  in  its  simplest  mode. 

4.  The  three-dimensional  array  was  reduced  tc  a  planar  horizontal  array.  Here 
we  are  returning  to  r‘-e  original  array  st rue tui  e  —  the  two  dimensional  arrey. 

Su:;— ary  results  ‘rearing  out  expectations  are  shown  in  Table  3,  The  maximal 
ratio  combining  scheme  gives  an  overall  moan  SNR  of  at  least  6  dB  better  than  the 
other  arrangements  except  for  EGC  case  which  is  a  close  second.  Interestingly, 
the  ratio  of  c (SNR) /<SN">  is  much  smaller  in  the  diversity  modes  than  in  the  two- 
dimensional  case.  This  too  is  expected.  The  diversity  inodes  arc  partly  effective 
in  resolving,  the  multipath  and  e.voidirg  the  nou-cohercnt  interference  of  the  multi- 
path  cr.'io.ii  ills,  In  the  two- d  inens Jo;,  ’ L  case  all  rays  entering  the  relatively  wide 
vertical  array  be,. , -.width  are  combined  non-cohert  nt  ly.  There  is,  therefore,  consid¬ 
erable  amplitude  ••aviation  dy ■  nd lug  on  the  relative  phases  of  the  accepted  rays. 

In  the  fixed  hn,.;;i  thre  -  dm  .  nsfon.il  rare  there  is  also  a  high  ratio  c  (SNR) /<  SNR> 
presumably  r.  re.-u’.t  of  the  narrow  vertical  hear  width  which  may  or  may  not  see  ar¬ 
riving  acoustic  energy. 

An  important  property  oi  the  array  will  be  its  response  to  sources  off  the 
azimuth  of  focus.  Tf  we  wore  i<>  inng.im  swing,  iv,  the  focus  away  from  a  source  which 
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MRC 

EGC 

SC 

90°  Sector 

2--Dimensional 

<SNR> 

480 

367 

105 

32 

117 

a (SNR) 

182 

156 

48 

30 

107 

o  =  50^  for  all  cases 

y 

h  =  251  for  3-dimensional  arrays 

TABLE  3  COMPARISON  OF  SIMULATION  RESULTS 


neatly  fed  M  branches  independently  the  power  in  each  branch  would  on  average  drop  by 

a  factor  N,  the  number  of  elements.  But  the  final  output  with  MRC  or  EGC  being  the 

M 

coherent  sum  of  the  M  weighted  branches,  would  only  be  reduced  by  a  factor  —  from 
the  main  beam  power.  A  situation  of  this  sort  was  simulated  to  check  this  surmise. 
The  result  obtained  using  the  MRC  system  was  <SN’R>  =  182,  the  average  being  over  the 
same  set  of  random  variables  as  before.  This  figure  is  somewhat  below  (M/N)  (main- 
bean  <SNR>)  but  it  .indicates  that  these  techniques  do  exact  a  price  in  sidolote 
response.  In  this  calculation  there  was  no  signal  on  the  mainbenm . 

Further  study  of  si  delobe  effects,  with  and  without  a  main  beam  signal  present 
and  using  the  different  combining  schemes,  is  viewed  as  a  useful  next  step.  in 
addition,  methods  based  on  estimation  theoretic  principles  (e;g.,  maximum  likeli¬ 
hood  and  maximum  entropy  estimation)  should  be  considered  for  application  here. 

These  methods  inherently  maximize  on-target  sign;]  response  relative  to  off-target 
signals.  Applied  to  the  separate  beans  as  found  here,  or  even  to  the  entire  array, 
superior  sidelobe  rejection  characteristics  can  be.  expected. 


Fred  Haber 
Paul  Yeh 
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